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It  is  argued  that  Eulerian  foraulations  are  intrinsically  unsuitcd 
for  deriving  the  Kolaogorov  theory  because  low-order  Eulerian  aoBents  do 
not  express  sufficiently  well  a  statistical  dependence  of  noosiaultanaoas 
asiplitudes  that  accoapanies  the  convection  of  snail  spatial  scales  by  large 
spatial  scales.  Illustration  is  nade  by  applying  the  direct-interaction 
approximation  and  a  related,  higher  Eulerian  approxiaation  to  an  idealised 
convection  problem  and  to  a  aodified  Navier-Stokes  equation.  Convection 
effects  of  low  wavenumbers  on  high  wavenumbers  are  removed  in  the  aodified 
equation,  and  as  a  consequence  the  direct-interaction  approxiaation  for  it 
yields  the  Kolmogorov  spectrum.  Low-order  Lagrangian  moewnts  provide  a 
promisingly  more  complete  description  of  the  convection  of  small  spatial 
scales  by  large,  and  a  search  for  satisfactory  Lagrangian  closure  approxi¬ 
mations  seems  highly  desirable. 
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1.  INTROSCICriON  AMD  SUMMARY 

Th®  si®il«rity  th®ory  proposed  by  Koliwgorov^*^  for  the  tMll-scele 
structure  of  turbulence  is  intuitively  ettrective,  end  the  theory  has  received 
strong  experiwntal  support.  Yet  Kol^jgoroVs  ideas  have  resisted  deduction 
fro«  the  Navier-Stokes  equation.  The  present  paper  suggests  that  Culerian 
analytical  formulations  are  intrinsically  unsuited  to  the  task  of  deriving 
the  KolMgorov  theory,  and  it  ends  uith  a  plea  for  attacks  by  Ugrangian 
■ethods. 

The  argument  can  be  suBMrised  siaq>ly:  An  underlying  assvaption  of 
the  KolMgorov  theory  is  that  very  large  spatial  scales  of  motion  conv^t 
very  small  scales  without  directly  causing  significant  internal  distortion 
of  the  small  scales.  The  assuaption  usually  is  considered  to  be  consistent 
uith,  and  to  imply,  statistical  independence  of  small  and  large  scales.  3ut 
this  is  true  only  for  the  simultaneous  velocity  disf-ibution.  It  is  pointed 
out  below  that  simple  convection  of  small  scales  by  large  scales  i^lies 
statistical  dependence  in  the  many-time  distribution  and  that  this  dependence 
cannot  adequately  be  specified  by  low-order  Eulerian  moewnts. 

If  only  a  fm  Eulerian  moments  are  specified,  it  is  not  possible  to 
distinguish  between  two  quite  different  physical  situations.  In  one,  the 
Eulerian  time-dependence  of  the  small-scale  strootures  is  due  to  their  being 
swept  by,  undistorted,  by  the  large-scale  motion.  In  the  other,  the  time- 
dependence  is  due  to  intrinsic,  internal  distortion  of  the  small  scales. 

The  implication  is  that  no  closure  approximation  which  retains  only  low-order 
Eulerian  moments  uses  enough  information  about  the  flow  to  sort  out  convection 
effects  from  intrinsic  distortion  effects.  In  contrast,  low-order  Lagrangian 
moments  do  appear  to  provide  this  information. 
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TImm  eoasidcNtions  atom  out  of  an  invoitigatioo  of  th«  inartial* 
mi  disalpatioD'ranca  prodlctiona  of  tba  dlract-intaractioo  apppoaiaatioo^ 

•  falatad,  higbar  approalMtiea.**  A  nuabar  of  rosulta  froa  this  invatti* 
gatloe  ara  givan  la  tha  praaant  papar  in  ordar  to  build  and  support  tha 
argwant.^  Tha  diraet-intaraction  approaination  givas  a  salf-consiatant 
elaaura  of  tha  Cuiarian  statistical  aquations  at  tha  lowast  nontririal  laval. 
It  aaprassaa  triple  nonanta  in  taraa  of  covariancas  and  yields  cloaad  aqua* 
tioaa  for  tha  eoaarianoas.  It  ia  shom  balow  that  this  approxination  is 
iaoapabls  of  daseribing  adaquataly  tha  coovaetioo  of  snail  scales  by  large 
acalaa.  Tba  higbar  i^proaiaatioo  giaas  cloaura  at  tha  next  laval.  It 
avpraaaaa  higbar  noaenta  in  tarM  of  triple  noaants  and  yields  closed  aqua* 
tiana  for  tha  triple  nonanta.  This  approxinatioo  is  found  to  give  an  in- 
prevnd  but  still  aarioualy  iaparfact  rsprasantation  of  cooraction  effects, 
■aoauaa  of  tha  difficulty  with  coovactioo  affects,  naitbar  approKinatioo  is 
capabla,  without  nodificatico,  of  giving  correct  predictions  of  tha  inertial* 
range  spactnas  law. 

Thasa  nattara  arc  clarified  by  considering  a  nodifiad  Naviar*Stokaa 
aquation  ia  which  tba  ccnvacticn  of  any  given  spatial  seals  by  nuch  larger 
soaiaa  is  eonsistantly  ranoved.  Tha  nodificatico  nay  ba  intarpratad  as  an 
effective  transfomatioo  to  tba  quasi*Lagrangiaa  coordinate  systan  invoked 
by  Kolnogcrov.^  When  appUad  to  tha  nodifiad  aquation,  tha  direct-interaction 
apprenination  raprodueas  tha  Kolnogorov  scaUng  laws  and  inert ial-ranga  power 
law.  Tha  higbar  approOnation  offers  the  hope  of  yielding  a  reasonably 
•***“'*ta  ntaerical  prediction  of  the  Kolnogorov  inertial-range  constant  and 
tha  universal  Kolnogorov  dissipation-range  spectrun.  However,  this  asstnns 
in  advance  that  tba  Kolnogorov  hypotheses  are  valid.  It  does  not  justify  than. 


The  paper  coocludee  with  a  diacutsioo  of  a  Ugrangian  frMework  which 
seeu  capable  of  diacriainating  between  connection  effects  and  internal 
distortion  effects  even  at  the  covariance  level  of  statistical  description. 
The  associated  equations  are  eore  coaq>licated  than  the  Navier-Stokes 
equation,  however,  and  it  is  not  clear  that  si^>le  but  consistent  closure 
approxieations  can  be  fomulated. 


2,  STATCNCNT  OF  THE  KOLMOGOROV  HYPOTHESES 

The  basic  assumption  of  the  KolK>gorov  theory  Is  that  dyn«ical  effects 
significant  for  energy  transfer  do  not  exist  over  indefinitely  large  ratios 
or  spatial  scales  of  notion.  An  equivalent  statement  is  that  suck  dynamical 
effects  do  not  extend  over  indefinitely  large  waventaber  ratios  when  the 
velocity  field  is  Fourier-analysed.  The  physical  picture  is  that  large-scale 
motions  should  carry  small  eddies  about  -ithout  distorting  then.  It  is  not 
obvious  that  this  need  be  true,  but  the  idea  is  certainly  intuitively  pUusible. 

The  basic  asswptlon  suggests  the  following  more  detailed  hypotheses 
(but  probably  does  not  fully  imply  then  in  a  rigorous  way);  1)  The  transfer 
of  energy  from  small  to  large  wavenumbers  at  high  Reynolds  niabers  proceeds 
by  a  cascade  process  which  is  effectively  local  in  wavenumber.  2)  Detailed 
information  about  the  low-wavenumber  structure  is  degraded  and  lost  during 
the  cascade.  3)  In  consequence,  there  is  an  asymptotic  high-wavenia^r  range 
whose  structure  is  universal  and  depends  only  on  the  kinematic  viscosity  v 
and  the  rate  c  at  which  energy-per-unit-mass  is  cascaded.  H)  The  turbulence 
is  isotropic  in  this  universal  range,  regardless  of  anisotropy  at  low  wave- 


oiad>ers. 
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It  follows  that  tha  universal  energy  spectrum  has  the  for»^ 

m)  .  (2.i) 

Here  C  is  som  universal  function  and  k^  is  a  nominal  dissipation  wavenuiBber 
A  final  assi>ptioa  is  that  5)  there  exists  an  inertial  subrange  k«  k  in 
which  direct  effects  of  viscous  dissipation  are  negligible  so  that  E(k) 
depends  only  on  c.  The  inertial-range  spectrun  is  then 

E(k)  .  (2.2) 

An  extensive  discussion  of  the  Koliwgorov  theory  ai.u  its  iisplications 
is  given  in  Ref.  2. 


3,  AN  IDEALIZED  CONVECTION  PROBLEN 

It  is  not  obvious  that  Kolmogorov’s  hypothesis  about  convection  is 
correct  for  actual  flows.  The  ssall-scale  (high  wavenumber)  structure  of 
high-Reynolds  number  turbulence  appears  visually  to  include  narrow  but  quite 
long  vortex  filaments  and  sheets.  It  is  not  wholly  clear  that  structures  of 
this  kind  are  carried  about  by  the  large-scale  motion  without  internal  dis¬ 
tortion.  In  the  present  Section,  a  very  simple,  in  fact  physically  trivial, 
problem  is  posed  in  which  the  large-scale  structures  unquestionably  convect 
the  small  scales  without  distorting  them.  In  Secs.  4  and  5.  it  will  be  seen 

that  this  idealised  problem  presents  severe  difficulties  for  Eulerian  closure 
approximations . 

Consider  a  total  velocity  field  u  t  v  with  the  following  properties:  v 
is  constant  in  space  and  time  and  has  a  Gaussian  and  isotropic  distribution 
over  an  ensemble  of  realisations;  u  varies  in  space  and  is  very  weak  compared 
to  vj  at  time  t  «  0,  u  is  statistically  independent  of  v  and  has  a  Gaussian, 


homogeneous,  isotropic  ensemble  distribution.  Assume  that  viscous  effects 
and  terms  bilinear  in  u  both  may  be  neglected  in  the  Navier-Stokes  equation. 
Then  any  Fourier  component  of  u  satisfies 


3u(k.t)/3t  =  -i(k‘v)u(k.t),  O 

whence 

u(k.t)  =  e'"!^*^^u(k,0).  (3.2J 

All  moments  of  the  u  field  can  be  obtained  from  (3.2)  and  the  stated 
statistical  properties  of  v  and  u(k.O).  For  example,  the  time-correlation 

R(kjt.f)  =  <u(k.t)«%.f)>/[<|u(k,t)|2X|u(k,f)|^i/2  (3.3) 

has  the  value 

where  v^  is  the  rms  value  of  v  along  any  axis.  The  last  member  of  (3.4)  is 
obtained  by  expanding  the  exponential  in  the  average,  using  the  rules  for 
evaluating  moments  of  a  Gaussian  distribution,  and  then  resumedng  the  series. 
Any  other  desired  moment  can  be  obtained  similarly. 

At  t  =  0,  the  u  and  V  fields  are  statistically  independent  by  assumption. 
Since  v  gives  simple  translation  of  u,  clearly  the  simultaneous  values  of 
the  two  fields  are  also  statistically  independent  at  any  later  time.  However, 
the  many-time  distributions  are  not  independent.  The  joint  distribution  has 
nonzero  cumulants  of  all  orders.  For  example,  it  follows  from  (3.2)  that 

^v.u.(k.t)u;(k.f)>  =  .ivjk.(t.f)exp[.  i  v3\^t-f)^Xu^U.0)u;(k.0)> 

(3.5) 

Similarly,  it  may  be  verified  that  the  more  general  c(*ulant  of  the  fora 

^v.v^...v^u.(k.t)u*(k.f)>.<v.v^--v,Xu.(k.t)u*(k.f)> 


(3.6) 
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dMS  not  vanish  unlaas  t  •  t*.  Thasa  rasults  illustrata  an  i^tortant  fact: 
Slapla  convactioD  of  saail-aeala  flea  coi^>OQants  by  larga-scaia  eon^>onants 
doaa  not  inply  statistical  indapandanca  of  tha  two  scalas.  On  tha  contrary, 
it  raquiras  a  statistical  dapandanea  which  involves  eusnilants  of  all  orders. 

It  has  bean  assuMd  so  far  that  tha  initial  u  field  is  purely  Gaussian, 
which  iaplias  statistical  indapandanca  of  different  Fourier  w^litudas.  Now. 
sui^osa  that  at  t  «  0  an  incraaant 

d«j(k.O)  •  i*'^<«jn-k^k„A^)u*(p,0)u*(q,0)  (3.7) 

is  added  to  tha  aapUtuda  for  a  particular  k,  whore  p  and  q  are  two  wave- 
vectors  such  that  k^£tq«0.  As  before,  let  u(p,0)  and  u(q,0)  be  statis¬ 
tically  indapandant  of  aich  other.  Tha  factor  ik  -k.k  ensures  that 

■  j  n 

Au(k,0)  is  incosfirassibla  and  real. 

Consider  tha  triple  correlation 


S<k,p„q;t,t',t") 


(kiUi(p,t»)u^(q.t")Au^(k,t)> 

<^kiUi(p,0)Uj(q,0)AUj(k,0)> 


(3.8) 


Tha  quantity  Au(k,t)  ovolvas  accordin*  to  (3.2).  It  than  follows  readily 
that 


S(k,p,q;t,t',t")  »  <a' 


(3.9) 

This  triple  correlation  is  of  a  kind  which,  in  actual  turbulence,  would  be 
associated  with  energy  transfer  «oog  tha  sodas  k.  p.  and  q.  Note  that  on 
tha  diagonal  t  »  t*  .  t"  tha  correlation  is  indapandant  of  t  [since 
k  ♦  p  ♦  q  .  0].  This  is  a  forsal  expression  of  tha  fact  that  translation 
by  tha  unifom  v  field  does  not  distort  tha  u  field.  Away  fros  tha  diagonal, 
tha  translation  does  induce  a  tine -dapandanea,  as  it  does  for  tha  double 
correlation  R(kit,t’). 


I  Th«  for  R(k;t,t')  «nd  S(k,p,q;t.t' .t")  hav«  bean  obtained  under 

I  the  aaaunptioo  that  ▼  ia  strictly  constant  in  space.  Ho*»erer,  these  results 
i  also  represent  asy^totic  lisdts  for  a  nore  general  situation,  in  which  v  is 
confined  to  vary  ssall  (but  nonsero)  wavenuabers,  has  a  hoMgeneous,  isotropic, 
Gaussian,  statistical  distribution,  and  is  statistically  independent  of  the 
initial  u  field.  This  correspondence  is  fairly  clear  physically  but  not 
i,»^i*tely  obvious  natheMtically.  If  »  varies  slowly  in  space,  (3.1)  nust 
be  replaced  by 

au(k,t)/3t  «  -il^,k*v(k*)i!(k-l5',t),  (3.10) 

which,  in  contrast  to  (3.1),  couples  different  wavevectors  of  the  u  field. 
However,  it  is  not  difficult  to  verify  by  direct  calculation  (expansion  in 
power  series,  averaging  using  honogenelty  properties,  and  suning)  that  (3.9) 
and  (3.H)  still  are  true,  provided  that  the  following  conditions  are  satisfied: 
1)  The  quantity  <|u(k,0)l^>  is  a  snooth  function  of  k.  2)  Equation  (3.7)  is 
replaced  by 

Au  (k,0)  -  ik  '  \  ^ 

where  k  is  any  wavevector  within  sosm  finite  rolmm  of  wavevector  space  and 
the  station  is  over  all  p  and  q  which  fall  within  sow  other  finite  volvws 
of  wavevector  space  (the  three  voluws  not  overlapping).  3)  The  wavevectors 
excited  In  the  v  field  are  very  snail  coepared  to  k,  p,  and  q  and  also  very 
snail  coeq>ared  to  the  width  of  the  voluaws  invoked  in  condition  2).  4)  Tiws 

large  conpared  to  (v^k)"^  are  not  conridered.  Conditions  1)  and  2)  (dense 
distribution  in  wave-vector)  ensure  ergodic  properties.  Conditions  3)  and 
•*)  ensure  that  the  shear  associated  with  the  v  field  produces  negligible 


dietortion  of  the  u  field  in  x  space. 


-8- 


Thls  aor*  gsMral  tnaaplc  still  differs  in  two  ways  froa  real  turbulence. 
In  an  actual  flow  the  clean  split  into  low  wavenuaber  v  field  and  high  wav«- 
nuaber  u  field  is  not  possible.  Moreover,  the  assuaption  of  statistical 
independence  of  low  and  high  wavenuabers  at  toae  given  tine  is  not  clwarly 
justified  in  actual  flows. 


a.  DIR£CT-lhTCRACT10M  APPROXINATIOM 

The  direct'interactioo  approxiaation  has  been  described  in  a  series  of 
3  4  6*9 

papers.  '  *  The  basic  idea  aay  be  aost  sinplv  described  as  follows,  for 
the  case  of  hoaogeneous  turbulence:  The  energy  transfer  aeong  the  fourier 
aodes  is  associated  with  triple  correlations  aaong  triads  of  interacting  inodes. 
In  the  direct-interaction  approxiaation,  it  is  assvssed  that  the  correlation 
of  any  given  triad  of  Fourier  nodes  is  induced  by  the  continuous  direct 
dynaaical  interaction  of  the  triad,  acting  against  a  relaxation  process  which 
destroys  the  correlation.  The  relaxation  process  involves  two  contributions: 
viscous  decay  and  dynaaical  relaxation  due  to  interaction  of  each  of  the 
three  aodes  with  all  the  rest  of  the  Fourier  nodes.  A  consistent  analv^.'ical 
expression  of  these  ideas  yields  fomulas  for  triple  correlations  in  terms  of 
two  functions,  the  velocity  covariance  and  the  average  response  function  of 
a  Fourier  aaplitude  to  inf initesiaal  disturbances.  The  approximation  then 
yields  closed  integrodifferential  equations  which  determine  these  functions. 

Niaerical  integration  of  the  direct -interact ion  equations  has  suggested 
that  the  approxiaation  gives  a  fairly  adequate  quantitative  description  of 
the  decay  of  isotropic  turbulence  at  modest  Reynolds  numbers.^  At  high 


Reynolds  numbers,  the  approxiaation  yields  a  universal  high-k  coectrun 


9 


E(k)  » 

0  o 


(k.l) 


where  v  is  the  r«s  value  of  velocity  in  any  direction  and  C*  is  a  universal 
0 

function.  The  dissipation  wavenuaber  k^  is  related  to  the  Kolaof^orov  wave- 

nuaber  k  bv 
s 

kd  *  C*.?) 

2 

where  is  the  Reynolds  nuaber  based  on  v^  and  the  Taylor  aicroscale. 

The  inert ial>range  spectrua  corresponding  to  (s.l)  is 

E(k)  *  (cvQ)^^^k*^^^C'(0).  (k.3) 


It  was  deduced  in  Ref.  3  that  the  direct -interact ion  anoroKiaation 
yields  a  local  cascade  of  energy,  as  called  for  by  Kolaocorov*s  assiaotioos. 
by  this  it  is  aeant  that  contributions  to  the  energy- transfer  into  an  inertial 
^  wavenia^r  k  froa  nodes  p,  q  is  negligible  if  the  ratio  of  any  pair  of  the 

wavenvasbers  k,  p,  q  is  large.  The  discrepancy  between  (4.1)  and  (2.1)  arises 
froa  a  rather  subtle  effect.  According  to  kclaogorov,  the  excitation  in  the 
energy-containing  range  does  not  affect  the  internal  dynaaics  of  the  energy- 
transfer  process  at  high  wavenuabers.  In  the  direct -interact ion  approxinat ion , 
the  rate  at  which  local  transfer  takes  place  at  high  k  depends  on  the  nagnitude 
v^  of  the  excitation  in  the  energy  range.  This  is  because  the  energy-range 
f  excitation  turns  out  to  contribute  to  and,  in  fact,  doainate  the  relaxation 
I  process  which  destroys  triple  correlations  aaoog  high  wavenuabers. 

Recent  experiaental  evidence^^  gives  strong  supoort  to  (2.1)  and  rules 
out  (4.1)  as  a  correct  asyaptotic  law.  In  the  present  Section,  theoretical 
support  for  this  conclusion  is  obtained  by  aoplying  direct -interact ion  to 
the  idealised  problem  of  Sec.  3.  The  approximation  is  found  to  yield  a 
relaxation  effect  of  the  v  field  on  triole  correlations  of  the  u  field  which 
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h«r«  is  clsarly  spurious.  This  doss  not  sstsblish  ths  Kolsogorov  thsory, 
but  it  doss  sesn  thst  ths  rssult  («.l)  should  not  bs  considsred  svidsncs 
ajcsinst  ths  thsory. 

In  ordsr  to  fonsulats  ths  soproxinstiofi  for  ths  idsalizsd  convsction 
problsB,  it  is  Dscsssary  to  introducs  ths  assra^s  rssoonss  function  G(k;t,t') 
of  ths  Tourisr  sods  k  to  infinitessal  psrturbations.  This  say  bs  defined , 
in  ths  prsssnt  case,  by 

«^^G(kjt,t')  *  <Cj^(k;t,t' )>, 

Cjj(k;t,t')  •  6u^(^t)/«f^(k,t'),  Cj^(k;t,t')  «  0  (t«t*).  (•».•*) 

whsrs  f  is  a  solsnoidal  forcing  tern  introducsd  on  ths  rirht-hand-sids  of 
0.1)  and  A/6  deiiorss  functional  diffsrsntiation.  (Sss  ^fs.  3  and  7.) 

Ones  ioducsd,  6u  is  sinply  translatsd  by  v.  Ksnes  ths  sxact  valus  of  G  is 

G(kit,f)  .  R(k;t,f)  (t  >  f).  (-.S) 

Ths  dirsct'intsraction  aporoxinatians  for  G  and  R  have  bssn  obtained 
orsviously  for  pmblsas  formally  identical  with  ths  prsssnt  one.  Ths  results 

ars 

G(k;t,t')  «  J^t?Wjjk(t-t')J/VQk(t-t')  (t  >  f).  (•».&) 

Rikjt.t*)  «  C(k;t,t*)  (t  >  t’),  (“.T) 

where  is  ths  first»ordsr  Bssssl  function.  This  rssult  for  G  can  bs 
represented  as  a  particular  partial  sv«Mtian,  to  all  orders,  of  ths  series 
obtained  by  expaodin*  ths  middle  member  of  («.S)  in  oowsrs  of  t-t*  and  then 
swsranioK.^  iUjuation  (k.6)  can  be  derived  in  either  of  two  wavs.  If  the 
V  field  ia  strictly  constant  in  space,  as  ass^Msd  at  ths  bs^inninr  of  Sec.  3, 
ths  problem  bscoMS  focsMlly  identical  with  ths  random  oscillator  orobles 
treated  in  Raf.  k,  and  ths  direct-interaction  results  are  obtaixtsd  bv 
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considering  «  collection  of  systems.  If  the  v  field  contains  sieail  nonzero 
wavenvaihers,  the  sane  results  are  obtained  by  the  more  physical  analysis 
of  .^f.  3. 

The  direct-interaction  approxieatian  for  S( k ,p ,q;t ,t '  ,t" )  is 

SiJifP.qit.t' ft")  *  G(k;t  ,0)K(p;t  •  .0)h(q;t",0).  (“.B) 

This  result  also  is  obtained  for  either  choice  of  v  field.  For  present 
purposes  assume  that  the  v  field  contains  nonzero  wavenui^ers.  The  initial 
correlated  increiaent  du^(k,0)  ftiven  by  (3.11)  nav  be  considered  the  result 
of  a.1  iepulsive  interaction  anon^  the  u  field  modes  at  t  *  0.  For  t  >  0, 
there  is  no  further  direct  cynaeical  interaction  aeionc  k,  p,  and  q.  The 
:.rocedure  of  Ref.  3,  Sec.  2  yields 

<u.(p,t’)u.(q,t")iu.(k,t)>  «  ik  («  -k  k  /k^) 

1—  j—  erurn 

■  <^“:i?.t)u*<P.0)uJq.t")u*(q,0)C  (k;t,0)>  *  ik  (6  -k  k  /k^) 

i-  n-  -j  —  j*—  emrn 

-  <Ui(p.t)u*(?,0))<u^(^.t")u*(q,0))^C^^(k;t,0)),  (e.S) 

whence  {•*,&).  The  factorins  of  the  average  is  not  oart  of  the  aporoxieation 

,  3 

Put  is  an  exact  consequence  of  hosiof eneity .  The  direct-interaction  approxi- 

aation  is  invoked  only  in  asserting  the  first  equality  in  (•».9).  The  approxi- 
■ation  excludes  the  "indirect"  contriPutions  to  S(k ,p,q ;t ,t * ,t" )  which  arise 
iron  initial  increaents  ^j(k-k',0),  where  k*  is  a  wavenuaber  excited  in  the 
V  field.  If  no  approxieation  is  made,  the  interaction  with  v  induces  a 
statistical  dependence  of  du^(k,t)  on  Au^(k-^*,0)  and  a  dependence  of 
Ui(p,t*  )u^(q,t")  on  u^(^k' ,0)u^(q,0)  and  Uj(p,0)Uj(qek  '  ,0).  In  addition, 

■ore  coaplicated  dependencies  are  induced. 

To  what  extent  are  the  direct-interaction  results  faithful  approximations 
to  the  exact  3,  P.,  and  S  functions?  The  results  for  3  and  K  aay  be  considered 


12 


qualitatiwly  accaptabl*.  Tbay  fi*«  ralaxatioo  and  corralatiao  tinas  ''(v^k)*^, 
in  accord  irith  tba  auct  rasulta.  Tba  principal  flaa  in  tba  diract-intaractioo 
aapraaaioDS  la  snail  daapad  oaclllatioaa  sbicfa  haaa  no  counterpart  In  tba 
aaact  aspraasiona.  (Saa  Tig.  2  of  Raf.  3.)  Tba  appraalaatioo  for  S  showa 
a  such  nora  sarlous  fault.  Tba  aaact  function  is  constant  on  tba  diagonal 
t  *  t*  ■  t",  utalla  tba  dlract'intaraction  result  decays  on  tba  diagonal. 

Away  fron  tba  diagonal,  tba  qualltativa  babaaior  of  tba  two  expressions  is 
slnilar. 

For  t  «  t*  ■  t",  S  is  tbs  kind  of  triple  corralation  which,  in  actual 
turbulence,  is  associated  with  anargy-transfar.  It  was  pointed  out  earlier 
in  this  Section  that  whan  diract-iataraction  is  applied  to  actual  turbulence, 
ralaxatioB  affects  aasociatad  with  anargy-ranga  convection  account  for  the 
diffaranca  between  (2.1)  and  (b.l).  Tbasa  affects  are  precisely  analogous 
to  the  spurious  decay  of  S(k,p,q;t,t,t}  glean  by  (k.8)  in  tba  present  axaapla. 
Tba  exact  natura  of  tba  iataraction  between  bigb  and  low  waeanunbars  in 
actual  bigb  Reynolds  suiter  flows  cannot  be  inferred  fron  tba  present  con- 
aidaratioQS.  Howsear,  it  would  be  niraculous  if  tba  spurious  affects 
associated  with  diract-interaction  bare  ware  sonabow  exactly  coapansatad  in 
actual  flows.  Henoa  there  is  no  justification  for  considering  tba  result 
(b.l)  as  aeidanca  against  tba  Kolaogorov  bypotbasas. 


S.  A  tUGHCR  CULERIAN  APPROXINATIOR 

Tba  diract-intaractioo  approniaatioo  closes  tba  Cularisn  statistical 

aquations  at  tbs  lowest  nootrieial  laeal.  It  proeidas  a  fomula  for  expressing 

a  ^ 

triple  Bonaats  in  tams  of  tba  ealocity  cowarianoa  <u(k,t)-u  (k.t'pand  tba 
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awrac^  raspcns*  ftactloo  G(k;t,t*)  of  an  Isotropic  floa  field.  The  final 

equations  contain  only  the  latter  two  quantities.  Tor  the  prohlea  of  Sec.  3, 

this  cloaure  (ires  a  qualitatively  acceptable  approaiaatioo  to  the  do«k>le 

correlation  R(k;t,t')  but  a  seriously  faulty  representation  of  convection 

effects  CX-.  the  triple  correlation  S.  The  question  naturally  arises  of 

whether  this  situation  can  be  raaedie4  by  a  higher  closure  approx^sation  in 

which  the  final  equations  are  at  the  level  of  triple  nonents.  An  apparently 

self-consistent  approaiaatioo  of  this  kind  was  introduced  in  a  previous 
s 

paper.  When  fomulated  for  isotropic  turbulence,  the  final  equations 
involve  the  following  three  kinds  of  higher  statistical  quantities: 

(uj(k,t)u^(p,f)u^(q,t-)>.  <^u^(k,t)6Uj(p,f)/6f*(q,t-)>, 

\«^u^(k,t)/«fJ(p,t‘)«f*(q^,t-)N, 

where  k  ♦  p  f  q  «  0  and  f  is  an  arbitrary  infinitesiaal  forcing  tera  added 
on  the  right -hand-side  of  the  Mavier-Stokes  equation.  An  extended  treatnent 
of  the  higher  approxlaation  is  in  preparation.  Sane  results  obtained  for 
the  idealised  convection  problea  of  Sec.  3  will  be  stated  now  and  coapared 
with  the  direct -interact ion  results. 

The  exact  equation  for  G(k;t,t')  in  the  idealised  problea  is 

3G(k;t,t')/>t  »  H(k;t.t*)  (t  >  t’),  G(kit,t)  »  1,  (5.1) 


where 

B(k;t,t')  «-i<it.*v*u^(k,t)/«f.(k,t*)>  (not  siMd  on  j)  (5.2) 
if  V  is  conatant  in  space, or 

H(kit,t')  »  -i[^,s^-»(k’)4u^(k-k',t)/4f^(k,t')>  (5.3) 

if  V  contains  very  snail  nonzero  waveniabers.  Here  f  is  an  azoitrary  farcing 


t*VB  addad  to  tiM  ri^bt'lMod-sido  of  0.1)  or  0.10).  It  is  convonisnt  to 
latroduoo  tlM  Lip  loco  troMfora 

Mk.i)  »  /  s‘***‘***6(k;t.t')il(t-t').  (5.S) 

0 

TIm  tromforo  of  (S.l)  is 

iC(k,i)  ■  1  ♦  Nik,!).  (5.S) 

TIm  dlroct-iatsroctioB  cinuro  forouli  for  U(k,s)  is 

i(k,i)  .  -sjk^t«k,s)]^.  (5.6) 

If  (S.6)  is  subotitutod  into  (S.S),  tboro  rosults  tiM  tmsfoiM  of  (s.6). 

Tbs  kifbar  opprooiostion  yislds  tbs  cloours  squitioo 

H(k.i)  .  -(Sj,k)2[ak,i)f.  (o^,k)'^tG(k,s)^^tH(k,•))^  (5.7) 

•ots  tbit  (S.6)  sj^rsssss  H  os  i  fisactioo  of  G  sbils  (5.7)  cootiios  H  oo 
both  sidos  of  tbs  oquitioo.  Cquitions  (S.5)  sod  (S.7)  uoy  bo  solsod  tofstbor 
to  yield  G(k,i).  Tbs  rssult  bos  boon  obtiinsd  sad  discussed  prsviously. 

Tbs  hifbsr  opproKiastioo  for  R(kit,t')  sitisfiss  0.7),  ia  nrn— iT  with 
both  tbs  suet  sad  dirsct'iatsrsctioB  rssults. 

Lot  ^ 

S(k,p,^jl,b,c)  «  ///  s“^"***’"^"s(k,p,q;t,t*,t")dtdfdt".  (5.0) 
0  -  -  - 
A 

Tbs  diiuct'iatsractioa  closurs  fosMuli  for  S  is 

S(k.p,q:i,b,c)  >  G(k,s)R(p,b)R(q,c),  (5.9) 

tbs  trsBSforo  of  (i.S).  Tbs  bigbsr  opproiiaatioo  for  S  sitisfiss  tbs 
soMsbit  coaplieitsd  squitioa 

S(*[.P.q»**b,c)  »  G(k,s)R(p,b)R(q,c)  ♦  ^k,j>,qis,b,c)S(k,p,q;s,b,c) ,  (5.10) 


sbsrs 
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A  •  ^  A  ^ 

D<k.P.qi*.t>,c)  «  *»(p.b)T(q,c)*yy(q.c)T(k,«)»xT{k,«)T{p,b),  (5.11) 

T(k,«)  .  (»jjk)-^(k.«)/a(k,«)  .  (»ok)-^{[G(k,.)]'^.*\. 


{LG(k,*)J  -a).  (j 

Th.  quMtitie*  X,  y,  1  *r«  th*  cmIms  of  th*  intorior  angUs  opposit*  k, 

P.  q,  ra«pocti*«ly,  in  a  triangle  wfaoae  sidaa  are  k.  p,  q.  Equation  (5.9) 
yielda 

'  G(k,a)R(p,b)R(q,c) 

s(k.P.5i..b.c)  . 

1  -  D(k,p,q;a,b,c) 


It  will  be  noted  that  the  higher  approxinationa  for  H  and  S  consist  of 
the  direct-interaction  contributions  plus  teuM  which  involve  H  and  S  th  ■»- 
wives.  The  added  terM  correspond  to  the  inclusion  of  additional  infinite 
subclasses  of  terw  fr«»  the  expansions  of  the  exact  H(k;t,t')  «d 
S(k,p,q;t,t' ,t")  in  power*  of  the  tine-argi»ents.  In  th*  diagrw 
of  Raf.  the  added  terw  -onstitute  iterated  vertex  corrections. 

To  what  extent  is  the  higher  approxination  an  iaprovenent?  The  naw 
values  for  G  and  R  turn  out  to  be  very  close  approxiMticns  to  the  exact 
values  and  represent  a  narked  quantitative  inprovewnt  over  the  direct- 
interaction  approxination  (see  Ref.  a,  rig.  13).  The  question  of  prin. 
interest  here  is  the  behavior  of  S  on  the  diagonal  t  =  t*  =  t",  which  na-r 
he  explored  by  taking  the  inverse  transfer*  of  (5.13).  The  analysis  is 
ciabersone  and  leads  to  n«rical  evaluations,  tiowever.  an  estiaate  of  the 
i^rovewnt  in  the  behavior  on  the  diagonal  nay  be  obtained  easily.  In  all 
three  cases,  exact,  direct -interact ion,  and  higher  approxination, 
S(k.p,q;t,t’,t-)  decays  away  fro*  the  diagonal  in  nor*  or  less  the  sm*  way. 
Thu*  th*  behavior  on  the  diagonal  nay  be  crudely  estinated  fron  th*  value 
of  the  integral  S(k.p,q;0 .0 .0) .  For  the  exact  S.  this  quantity  is  infinite 
because  there  is  zero  decay  on  the  diagonal.  In  th*  direct -interact ion 
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HyroKi—tion, 

^;0,0,0)  *  1.  {5.1») 

A 

la  tlM  bicbar  approaiaatioe,  S(k,p,q;0,0,0)  rcMins  fiaitt,  which  show*  that 
dacay  on  tbs  still  occurs.  Howsvcr,  than*  is  on  increosa  orsr  tha 

diract-intsraction  volus,  which  io<licatM  that  th«  decoy  is  slower.  The 
factor  of  lacrease  dspeads  oa  the  shop*  of  the  triongls  foraed  by  k,  p,  and 
4- 

T^la  I  shows  the  eslusi.  of  S(k,p,q;0,0,0)  given  by  (S.13)  for  several 

A 

choioes  of  waveaiartMr  ratios.  The  increase  of  S(k,p,^;0,0,0)  coaes  frca 
two  factors,  first,  G(k.O)  Is  bigger  in  the  higher  approxiastion  [the 

A 

direct -late  recti  on,  higbei^approitiaatian,  and  exact  values  of  v^kG(k,0)  ore, 
respectively  l,Ul  ♦  «  1.272,  and  (1/2)^^^  »  1.253J.  The  second 

factor  Is  the  denoalaator  in  (S.13).  This  factor  olons  is  actually  a  acre 
valid  asasure  of  tbs  dscay  of  S(k,p,q;t,t,t)  than  is  S(k,p,q;0,0,0)  itself. 

The  lat_ter  quantity  reflects,  in  part,  changes  in  rete-of -decoy  perpendicular 
to  the  diagonal. 


TAhLC  I.  Values  of  S(k,p,^;0,0,0)  in  the  higher  opproxiaatian. 


k:p;q 

1/Cl-D<k,p,q;0,0,0)) 

v^pqS(k,p,q;0,0,0) 

1:1:1 

13.71 

2B.21 

1:2:2 

6.66 

13.70 

l:s:s 

h.OS 

8.3a 

2.62 

S.39 

3:S:S 

7.e2 

15.27 

1:1:^ 

7.9a 

16.3a 

doother  neoeurs  of  the  li^>roeaaeat  ia  the  higher  approxiastion  is  the 
behavior  of  S(k,p,qit,t,t)  at  t  a  0.  In  the  direct -interact ion  opproxiaation 


Mcood  <l.riY«ti»«  with  respect  to  t  is  negstiee.  In  the  higher  spproxi- 
■stioo,  it  is  not  hsrd  to  show  thst  the  first  three  derivstires  eanish  end 
the  fourth  derivative  is  negative. 

In  the  present  context,  the  fact  of  principal  interest  is  sieply  that 
the  spurious  reUxatioo  effects  of  the  direct-interaction  approxination  are 
reduced,  but  not  elininated.  Direct-interaction  provides  closure  at  the 
covariance  level  and  gives  a  qualitatively  satisfactory  double  correlation 
R(k;t,t').  It  night  be  thought.,  by  analogy,  that  the  higher  apprcxination, 
which  gives  closure  at  the  triple  nosKr.t  level,  would  give  a  qualitatively 

satisfactory  triple  correlation  S(k,p,q;t,f ,t-).  But  clearly  this  is  not 

the  case. 


6.  .'WDiriED  BAVIER-STOKXS  EQUATION 

The  Navier-Stokes  equation  can  be  altered  so  as  to  renove  consistently 
the  convection  effects  of  low  wavenunbers  on  high  waven;n*ers.  The  nodi  find 
equation  serves  to  give  a  precise  statenent  of  Kolnogorov's  basic  hypothesis 
and  to  confim  the  role  of  convection  effects  in  producing  the  discrepancy 
between  the  direct-interaction  and  Aolnogorov  spectra, 
nodified  equation  nay  be  writven 


{>/»ttvk^)u(k,t)  *  -il*.[u(k-k').k']u(k*,t)-ikv(k,t),  k.u(k,t) 

where  »(k,t)  is  the  Fourier  anplitude  of  kinenatic  pressure  and 
tams  such  that  " 


=  0,  (6.1) 
onits  all 


|k-k' I  <  k/a  or  !k-k*l  <  k'/a. 

The  paraMter  a  is  an  arbitrary  cut-off  ratio  >1.  The  case  a 
unnodi f led  Navier-Stokes  equation. 
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rourUr  transfora  of  ths  convection  operator 

& 

u*V  .  To  see  the  significance  of  trwtcating  this  sua,  consider  a  given  wave- 
vector  triad  tkt  Jft  that  *  0  and  aithar  q  <  k/o  or  q  <  p/o. 

The  intaractioo  of  this  triad  is  r. 'presented  by  a  set  of  teras  in  the  Navier- 
Stokes  equations  for  u(k,t),  u(p,t),  and  u(q,t).  The  tenw  in  this  set  which 
are  onittad  and  retained  in  the  Modified  equations  are  shown  in  Table  II. 

The  oaitted  and  retainad  tarns  saparately  conserve  the  energy  of  t*-*  triad. 
Thus  the  nodified  equations  are  energetically  consistent.  The  teras  onitted 
correspond  to  the  convection  and  straining  of  snail  spatial  scales  by  large 
scales.  It  is  clear  fron  Tabla  II  that  thesa  tarns  are  associated  with 
energy  transfer  out  of  q  into  k  and  p»  since  the  teras  do  not  contribute  to 
the  of  u(q,t)  in  tne  uanodified  equations.  However,  these  teras  do 

reprasant  a  straining  action  which  can  yield  energy  transfer  between  and  p. 
The  teras  retained  represent  a  transport  of  the  nonentun  of  the  large-scale 
notions  by  convective  action  of  the  snail  scales  (snail  and  large  here  are 
relative  teras).  These  terns  conprise  the  eddy-viscosity  nechanisn  (v  altered 
fron  the  original  Navier-Stokes  equation)  by  which  energy  is  transfered  fron 
q  to  k  and  p. 


TA8LE  11.  Tenn  onitted  and  retained  in  the  nodified  Mavier-Stokes 
equation  when  q  <  k/«  or  q  <  p/*. 


Onitted 


In  k  Cq. 
In  p  £q. 
In  q  Cq. 


ilu( -^)-p]u(-p) 
i[u(-q).k  Ju(-k) 


Retained 


i[u(-p)‘q)u(-5) 

i[u(-k)-q]u(-q) 

i[u(  -k).p)\i(-p)ti[u(-p)  'kju^l-k) 


Hone 
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Kolaogorov's  basic  assuBption  is  that  tha  action  of  lar^s-scala  ■otions 
on  sufficiently  smaller-scalc  Motions  is  a  sinple  convection  process  that 
does  not  affect  appreciably  the  energy  dynaMics  of  the  snaller  scales.  This 
hypothesis  can  be  given  a  precise  statanent  in  terns  of  the  Modified  Navier- 
Stokes  equation:  If  the  cut-eff  ratio  a  is  taken  large  enough,  the  difference 
in  the  energy  spectrua  E(k)  for  the  Modified  and  unaodified  equations  can  be 
Made  as  snail  as  desired,  when  the  spectrua  is  noraalized  in  ter«s  of  the 
Kolaogorov  similarity  parameters.  For  given  a,  this  difference  is  independent 
of  as 

Since  the  modified  equation  removes  the  convection  effects  of  large 
scales  on  small  scales,  it  represents,  in  effect,  a  transformation  to  the 
quasi-Lagrangian  coordinates  called  for  in  the  original  statement  of 
Kolmogorov's  theory.^  In  this  regard,  the  modified  equation  has  the  ad¬ 
vantage  of  being  self-consistent  for  all  t.  It  is  difficult  to  construct 
s  well-defined  transformation  to  locally  co-noving  coordinates  in  x  space 
which  is  similarly  consistent. 

'he  direct-interaction  approximation  for  (6.1)  goes  through  just  as 
for  the  unmodified  Navier-Stokes  equation.  The  only  difference  in  the  final 
statistical  equations  is  that  the  geometrical  coefficients  which  occur  now 
have  altered  values.  Th.  new  values  are  given  in  the  Appendix.  When  a  is 
finite,  the  direct -interact ion  approximation  yields  an  assmptotic  spectrua 
of  the  Kolmogorov  form  (2.1)  rather  than  (4.1).  This  comes  about  as  follows: 
The  removal  of  cc'vection  effects  of  the  energy-ra ige  on  high  wavenumbers 
has  the  result  that  R(k;t,t')  and  G(k;t,t*>  for  large  k  have  decay  times 
determined  by  local  excitation  levels  rather  than  by  the  energy-range 
excitation  level  v^.  In  consequence,  the  relaxation  of  triple  correlations 
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moag  th.  high  i.  ind*p«nd*nt  of  v^.  The  Kolmogorov  spectrvj. 

then  resultt  inescepahly  because  it  turns  out  that  the  energy- transfer  is 
local,  as  in  the  us»odified  case.  R(k;t,t')  and  G(k;t,f)  may  nou  be 

interpreted  as  correlation  and  response  functions  measured  in  quasi-Lagrangian 
coordinates. 

Figure  1  shows  some  results  of  numerical  solution  of  the  direct- 
interaction  equations  for  several  values  of  a.  in  the  range  v  6C0-800. 

The  cv.rvs  were  obtained  as  follows.  Random  stirring  forces,  confined  to 
wavenumbers  below  those  shown,  were  invoked  in  order  to  give  a  statistically 
steady  state.  The  response  and  time-correlation  functions  were  fitted  to 
Gaussian  functions  of  difference  time,  thereby  yielding  closed  equations 
'Aich  contain  only  £(k)  and  the  characteristic  correlation  and  response 
times  R(k.O)  and  G(k.O).^  These  equations  were  discretized  in  1/’  octave 
k  bands  by  a  previously  tested  schesa’  and  solved  by  iteration.  The  overall 
numerical  error,  due  to  the  fitting  of  the  response  and  correlation  functions, 
the  discrstiting  in  k.  and  truncation  in  k,  is  estimated  as  less  than  10%  of 
local  value. 

The  function  plotted  is  the  one -dimensional  spectrum 

•  ^(k)  .  i  /^(l  .  k^/p^>p’Mp)dp,  (6,3) 

normaUsed  according  to  the  Kolmogorov  similarity  scaling.  Numerical  investi¬ 
gation  indicates  that  the  finite  o  curves  differ  inappreciably  frosi  the 
asymptotic  curves  for  in  the  range  of  k/k^  shown.  The  a  =  -  curve 

(unmodified  equation)  differs  inappreciably  from  the  =  ■  values,  if  it 
is  normalized  according  to  the  k^  scaling  associated  with  (4.1).  it  follows 
from  (4.2)  that  the  position  of  this  curve  on  the  present  Kolmogorov-scaled 
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Tuifl  the  high  wavenuabers  is  indspendant  of  v^.  The  Kolaogorov  spectrxa 
then  results  inescapably  because  it  turns  out  that  the  eneri^- transfer  is 
local,  as  in  the  unaodified  case.  R(k;t,t')  and  G(k;t,t')  may  now  be 
interpreted  as  correlation  and  response  functions  aeasured  in  quasi-Lagrangian 
coordinates. 

Figure  1  shows  soae  results  of  niaarical  solution  of  the  direct- 
interaction  equations  for  several  values  of  a,  in  the  range  %  600-800. 

The  curves  were  obtained  as  follows.  Randoa  stirring  forces,  confined  to 
wavenuabers  below  those  shown,  were  invoked  in  order  to  give  a  statistically 
steady  state. ^  The  response  and  tiae-correlation  functions  were  fitted  to 
Gaussian  f\nctians  of  difference  tine,  thereby  yielding  closed  equations 
which  contain  only  £(k)  and  the  characteristic  correlation  and  response 
tines  R(k,0)  and  G(k,0).^  These  equations  were  discretized  in  I/**  octave 

7 

k  bands  by  a  previously  tested  schem  and  solved  by  iteration.  The  overall 
nuaerical  error,  due  to  the  fitting  of  the  response  and  correlation  fxnctions, 
the  discretizing  in  k,  and  truncation  in  k,  is  estisksted  as  less  than  10%  of 
local  value. 

The  fraction  plotted  is  the  one-diaensional  spectr^a 

♦  .(k)  .  i  /  (1  -  k^/p^)p‘S:(p)dp,  (6.3) 

^  k 

noraalised  according  to  the  kolai^gorov  sieilarity  scaling.  Nuaerical  investi¬ 
gation  indicates  that  the  finite  a  curves  differ  inappreciably  froa  the 
asyaptotic  curves  for  R^^  *  •,  in  the  range  of  kA^  shown.  The  a  s  •  curve 
(unaodified  equation)  differs  inappreciably  froa  the  R^  *  *  values,  if  it 
is  noraalised  according  to  the  k^  scaling  associated  with  (4.1).  It  follows 
froa  (4.2)  that  the  position  of  this  curve  on  the  present  Kolaogorov-scaled 
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••tiMtts  of  C{0)  froM  diffartnt  sources  show  a  substantial  spread.  Grant, 
10 

Stawart,  and  Moilliet  have  obtained  C(0)  >"  l.**4  froa  aeasureaents  in  a 
tidal  channel  at  >  1000. 

The  higher  approximation  discussed  in  Sec.  S  does  not  remove  completely 

the  spurious  relaxation  of  high-wavenwber  triple  correlations  due  to  energy- 

range  convection  effects.  Conseqiiently,  it  leads,  like  the  direct-interaction 

-3/2 

Approximation,  to  an  asymptotic  k  inertial  range  law.  However,  the 
spurioua  effects  are  sufficiently  reduced  that  the  higher  approximation  may 
yield  a  rather  accurate  prediction  of  the  Kolmogorov  canstant  and  of  the 
Kolmogorov  dissipation  spectrum  if  it  is  applied  to  the  modified  Navier- 
Stokcs  equation.  For  this  to  be  so,  there  must  be  values  of  the  cut-off 
ratio  a  such  that:  1)  The  exact  normalised  E(k)  associated  with  modified 
and  unawdified  equations  differ  very  little.  This  requires  that  a  be  large 
enough  that  the  cut-off  removes  very  little  of  the  energy-cascade  of  the 
unmodified  equation;  it  is  presupposed  here  that  Kolmogorov's  hypotheses 
are  correct.  2)  The  errors  due  to  spurious  relaxation  are  sufficiently 
•■all  that  the  higher  approximation  gives  an  accurate  representation  of  the 
<iy&amics  of  the  eodified  eqmation;  this  requires  that  a  be  small  enough  that 
tba  raaiduai  spurious  effects  art  small  compared  to  the  real  relaxation 
effects  associated  with  energy-transfer  and  internal  distortion  among  local 
wavsniabera. 

An  investigation  now  in  progress  indicates  that  values  of  a  satisfying 
these  two  apposing  canditlons  should  exist  and  indicates  also  a  way  of  sati- 
mating  the  residual  errors,  on  the  aasimption  that  Kolmogorov  is  correct. 

The  difficulties  of  oumsrical  solution  poeed  by  the  final  aquations  are  severe. 
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7.  HC7TIVAT10H  FOR  LAGRAMGIAR  APPROACHES 

The  preceding  analysis  has  aade  clear  that  direct-interaction  gives 
spurious  relaxation  effects  by  low  eavenuaber  f lex^-coaponents  on  triple 
correlations  asong  high-uavenunber  coaf>ooeots.  It  is  pertinent  to  ask  vhy 
this  particular  kind  of  trouble  should  arise  in  an  approxination  which  is 
energetically  consistent  and  which  appears  to  provide  a  reasonably  faithful 
description  of  turbulence  at  Reynolds  nied>ers  low  enough  that  a  wide  range 
of  wavenusiers  is  not  strongly  excited.  In  the  author's  judgment ,  the 
underlying  trouble  is  elementary  and  can  be  expected  to  afflict  any  closure 
approximation  which  expresses  triple  correlations  as  a  function  of  E(k), 
R(k;t,t'),  and  G(k;t,t'). 

In  the  case  of  turbulent  interaction  among  wavenumbers  similar  in 
magnitude,  the  function  R(k;t,t')  exoresses  the  less  of  correlation  associated 
with  disordered  internal  distortion  of  the  flow  patterns.  The  niaerical 
success  at  low  Reynolds  numbers  suggests  that  this  internal  distortion 
produces  relaxation  of  triple  correlations  in  a  way  that  is  reasonably  well 
described  by  the  direct-interaction  closure  formula.  Dn  the  other  hand,  the 
energy-range  convection  at  high  Reynolds  numbers  produces  a  loss  of  correlation 
in  the  ai^>litude  of  a  high-waveni*ber  mode  which  is  not  associated  with 
significant  internal  distortion.  The  difficulty  here  is  that  it  is  impossible, 
given  only  R(k;t,t*),  to  say  how  much  of  the  correlation  loss  is  Cue  to  dis¬ 
tortion  and  how  much  to  simple  convection.  In  order  to  distinguish  the  two 
effects,  it  is  necessary  to  specify  the  higher  statistical  structure  of  the 
velocity  field.  In  particular,  it  is  necessary  to  say  something  about  the 
c\»ulants  of  form  (3.6)  which  describe  the  joint-distribution  of  low  and  high 
waveniabers.  This  higher  statistical  information  has  no  way  of  entering  the 
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dircct- interact  ion  scheaet  In  effect,  the  closure  forcula  maxes  tr.e  decision 
that  the  lose  of  correlation  in  individual  eode  aaplitudes  is  entirelv 
aseociated  with  intrir.sic,  internal  disordering  of  the  flow  sattems. 

Since  direct -interact ion  gives  qualitat ive ly  correct  behj.icr  for  the 
double  correlation  R(k;t,t')  even  when  convection  effects  are  strong,  it 
might  be  hoped  that  closure  of  the  Lulerian  statistical  equations  at  one 
step  higher  would  give  qualitatively  correct  triple  co-relations.  The 
results  prssented  in  Sec.  S  do  not  support  this  hope.  The  failure  eight  nave 
been  ecticipated.  In  order  to  specify  that  the  a  field  is  suffering  pure 
coovection  over  finite  times  of  evolution,  it  is  necessary  to  give  proper 
values  to  the  entire  infinite  set  of  ci*ulants  of  the  fore  (3.6),  It  is  net 
too  surprising  that  a  finita-order  truncation  of  the  statistical  equations 
fails  to  accomplish  this. 

Tba  analysis  in  this  paper  has  been  based  on  the  many-time  distribution 

of  the  velocity  amplitudes.  The  cumulants  (3.6)  describing  interdependence 

of  low  and  high  wavenwbers  all  can  vanish  if  the  amplitudes  are  measured 

sirultanaousiy.  iftiy  not  then  seek  a  closure  approximation  involving  only 

simultaneous  amplitudes  and  thereby  avoid  all  the  difficulty?  The  quasi- 

12 

oormelity  approximation  is  an  energetically  conservative  closure  of  this 
kind.  If  this  approxisMtion  is  applied  to  the  idealized  convection  probles 
of  Sec.  3,  it  is  easy  to  see  that  spurious  relaxation  of  triple  moments  of 
the  u  field  in  fact  does  not  occur.  The  trouble  is  that,  when  applied  to 
actual  turbulance,  the  approximation  also  gives  no  relaxation  of  triple 
correlations  produced  by  interaction  of  similar  spatial  scales.  Since  only 
simultaneous  amplitudes  enter  the  equations,  tne  approximation  has  no  way 
of  describing  loss  of  correlation  and  consequently  eskbodies  no  mechanise  at 
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ail  for  dynaaical  relaxatioo  of  dyiiaadcally  induced  correlations.  The  con¬ 
sequence  is  an  unph3rsical  overshoot  of  initially  built-up  enerfv-transfer 
and  an  eventual  prediction  of  negative  energies  in  the  spectrva  regions 
wnich  have  strongest  initial  excitation. Because  of  the  failure  to 
eabody  relaxation  effects  in  interactions  of  nearby  wavenuabers,  the  anproxi- 
aation  does  not  yield  Kolnogorov's  spectna. 

In  contrast  to  R(k;t,t'),  an  appropriately  constructed  Lagrangian 
double  correlation  appears  to  permit  a  clear  specification  of  whether  loss 
of  correlation  in  a  single  Fourier  node  is  due  to  internal  distortion  of 
the  flow  structures  or  to  siaple  convection  of  these  structsres.  L<et 
u(x,t|s)  be  defined  as  the  velocity  at  tiae  s  of  a  fluid  elenent  which 
arrives  at  the  point  x  at  tiae  t.  This  will  be  called  the  Lagrangian  velocity 
in  wnat  follows,  because  it  is  aeasured  along  the  trajectory  of  a  fluid 
particle.  If  t  is  fixed,  u(x,t|s)  nav  be  identified  with  the  Lagrangian 
velocity  as  sore  customarily  defined, with  x  serving  to  label  the  particular 
particle  whose  path  is  being  followed.  The  full  function  u<x,t|s)  may  be 
regarded  as  giving  an  infinite  set  of  Lagrangian  velocity  fields  (according 
to  the  restricted  definition)  with  each  tiae  t  specifying  a  different  labeling 
of  the  particle  trajectories. 

The  statistical  structure  of  u(x,t|s)  aay  be  specified  by  its  aoaents. 

If  the  field  has  zero  aean,  the  siaplest  nonvanishing  aoaent  is  the  covariance 


U.^(x,t;x' ,t' |».s' )  =  ^u^(x^t  |s)u^(x' ,t' Is' )')•  (7.1) 

This  function  has  four  tiae  arg\*ents  instead  of  the  two  which  appear  in  the 
Eulerian  covariance  ^u.(x,t )u^(x' ,t '  )^.  Consequently,  it  conveys  aore  infoe^ 
.  Consider  again  the  idealized  convection  orobler.  of  Sec.  3.  Under 


it  ion 
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th«  — wptiop  of  •potially  constant  ▼  field  and  negligible  terns  bilinear 
in  the  Lagrangian  ntlocity  is 

u(a,t|s)  3  u(s.t),  (7.2) 

since  the  uniforn  translation  by  jr  does  not  change  the  velocity  of  any 
particle.  Let  the  Lagrangian  tine-correlation  be  defined  by 

R(kit,f|».**)  *  <«(k.t|s).u*(k,f|*')>/l<|u(k,t|s)|^Xl!i<^f|s*)l^]^^^. 

(7.3) 

where  u(k,t|s)  is  the  Fourier  coefficient  of  u(ji,t|s).  By  (7.2), 

R(k;t,t'|s.s')  »  R(k;t,t’).  (7.**) 

According  to  (7.S),  R(k  ;t  ,t '  |  s  ,s*  )  exhibits  no  depexKience  on  s  and  s', 
the  tines  at  which  the  nelocity  field  is  neasured.  It  varies  only  with  the 
tines  t  and  t*  which  specify  the  particle  labeling.  On  the  other  hand,  if 
the  u  field  exhibited  internal  distortion  also  (change  with  tine  in  a  fraee 
noving  with  v)  there  clearly  would  be  a  dependence  on  t  and  t*.  Thus  it 
appears  possible  that  the  lowest-order  Lagrangian  tine-correlation  nay  convey 
enough  inforaation  to  decide  bow  nuch  of  the  change  in  given  Fourier  conponents 
of  the  velocity  field  is  due  to  sinple  convection  and  how  nuch  to  internal 
distortion.  This  suggests  that  it  is  worthwhile  to  search  for  a  consistent 
Lagrangian  closure  approxination  at  the  lewel  of  j(x,t ;x' ,t ' | s ,s' ). 

Uhfortiaiately ,  the  equation  obeyed  by  u(x^t|s)  is  nuch  nore  conplicated 
than  the  Hawier-Stokes  equation.  When  the  pressuxe  is  elininated,  the  Nawier- 
Stokes  equation  nay  be  written 

(3/Jt  ♦  u«7)u.(j^t)  «  wV^Uj(«,t)  ♦  ■i^(^)Cu^(*,t)u^(x,t)],  (7.5) 

wtterc 

I.  (V)  *-v"^J^/3x.3x  3x  , 
iSB  inn 

9*^w(x)  5  -(wv)*^/|,-K;|*^e<x')d^x' 


t 


(7.6) 


-27- 


foc  any  suitably  bahaved  function  »»(*).  The  field  u(x,t|s)  satisfies 


s 


(7.7) 


where 


t 


C(jt,t|r)  =  /  u(;,t|r')dr' 


(7.8) 


r 


■easures  particle  displacefwnt  in  the  interral  (r.t)  and  u(5,0|0)  is  the 


prescribed  initial  field.  Mote  that  u(*,t|T)  =  u(*,t)  by  definition.  The 
notation  e*?^  signifies  that  in  the  expansion  of  the  exponential  all  ^  factors 
appear  to  the  right  of  all  C  factors;  that  is,  ^  operate  on  The 

(Privation  and  interpretation  of  (7.7)  is  sieple.  The  operator  exp^C“]4*^3 
is  a  displacenent  operator;  e.g.. 


(7.9) 


exp^t-4*V)u(x,0)  s  if^**S*^^* 


Thus  (7.7)  states  that  the  velocity  field  at  later  tines  is  the  result  of 
self -convection  of  the  initial  velocity  field,  together  with  convection  of 
all  the  velocity  increnents  induced  at  later  tines  by  viscous  and  pressure 
forces. 

It  has  so  far  proved  difficult  to  carry  out  consistent  closures  of  the 
statistical  equations  generated  by  (7.7).  The  nonlinearity  is  transcendental, 
in  contrast  to  (7.5).  A  related  fact  is  that  the  inconpressibility  condition 
for  the  full  field  u(^t|s)  has  an  ankward  fore.  loccnpressibility  for  all 
tines  is  aut  optically  guaranteed  by  (7.7),  if  u(x,0|0)  U  divergence  less , 
but  it  is  difficult  to  nake  statistical  approxinations  which  preserve  this 
property. 

A  prelininary  application  of  the  u(x,tjs)  formalise  to  dispersion  of  a 

16 

scalar  field  by  honogeneous  turbulence  has  been  reported.  Closure  was 
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obtaiiMd  by  th«  crude  aoproximat ion  that  £(x,t|s)  is  ■ultivariate-Gaussiac. 

This  gives  SGBS  reasonable  predictions  for  the  evolution  of  the  scalar 
spectruB.  (The  results  are  quite  different  fro*  those  of  the  Eulerian  quasi¬ 
noreality  approxisLation. )  In  particular,  the  result  for  the  straining  effect 
of  large  eddies  on  snail  scalar-field  structures  seens  qualitatively  olausible. 
In  contrast,  the  direct-interaction  approximation  for  this  problec,^^  exhibits 
spurious  relaxation  effects  on  the  scalar  field  very  sinilar  to  those  dis¬ 
cussed  in  the  present  paper. 

The  assiaption  of  ■ultivariate-^:;aussian  i  is  exactly  valid  for  the 
idealised  problee  of  Sec.  3  ax>d  yields  both  (3.%)  and  (3.9).  The  reason  is 
that  vnder  the  conditions  iaposed  there 

Ux,t|s)  *  (t  -  s)v, 

and  v^  is  .Russian,  nouever,  this  closure  approximation  can  violate 
incoapressibility  and  thereby  lead  to  serious  troubles  if  it  is  applied  to 
(7.7)  vnder  eore  physically  interesting  conditions.  3ther  closure  aoproxi- 
mations  should  be  sought,  and  alternative  Lagrangian  formulations  should 
be  explored. 


ACM(OWL£3GO(CirrS 

I  an  grateful  to  I>r.  N.  n,  Gibson  for  supplying  tabulations  of  his 
original  data.  The  value  of  C*(0)  reported  in  Sec.  6  corroborates  a  calculation 
of  this  constant  made  several  years  ago  by  tr.  Kobert  <^roic)c,  using  a 
different  method. 

This  uorh  was  supported  by  the  Fluid  ^hmamics  Branch  of  the  Office  of 
Naval  Research  vader  Contract  Nonr  s307(00).  Facilities  for  the  numerical 
calculations  twrs  kindly  furnished  by  the  Goddard  Institute  for  Space  Studies. 


-29- 


APPEXDiX 

The  fio«l  isotropic  direct-interaction  equations  are  given  in  Sec.  3  of 
Ref.  3,  for  the  statistically  steady  case,  and  in  Sec.  2  of  Ref.  7,  for  free 
decay.  Wien  the  havier-Stokes  equation  is  nodified  to  renove  nonlocal  con¬ 
vection  effects,  the  geonetrical  coefficients  a(k,p,q)  and  b(k,p,q)  irtiich 
occur  in  the  final  equations  have  altered  values  if  any  pair  of  the  wavenunbers 
P»  q  *  ratio  which  exceeds  the  cut-off  paraneter  o.  There  are  no  other 
changes.  The  new  values  are 

«U,p,q)  *  a(k,q,p)  =  ^  (q^A^Xl  -  *^)(1  ♦  y^)  (k  >  oq  or  p  >  aq), 

a(k,p.q)  =  a(k,q,p)  =  y  (1  -  xy*  -  2y^«^)  (p  >  ak  or  q  >  ak), 

b(k,p,q)  =  y  (q^/k^Xxyx  ♦  xV)  (k  >  oq  or  p  >  oq), 

b(k,p,q)  =  2a(k,p,q)  -  b(k,q,p)  (k  >  oo  or  q  »  op), 

b(k,p,q)  *  a(k,p,q)  ♦  i  (x^  -  y^)  (p  >  ok  or  q  >  ok), 

where  x,  y,  and  z  are  the  cosines  defined  in  the  Refs,  and  the  oresent  text. 

The  elieination  of  v^  fro*  the  asyeptotic  inertial-range  functions  R(k;  t,t') 

2  2 

and  Glk;t,t*)  can  be  traced  to  the  factors  (q  A  )  in  the  expressions  above. 

There  are  also  other  ways  of  reswving  nonlocal  relaxation  effects  frcm 
the  final  statistical  equations.  Let  the  direct -interact ion  equations  be 
written  as  separate  equations  for  E(k),  R(k;t,t'),  and  G(k;t,t')  CLqs.  (0.2), 
(S.l),  and  (S.o)  of  Ref.  3j.  Let  the  E(k)  equation  be  unaltered,  but  in  the 
R  and  G  equations  eliainate  nonlocal  convection  effects.  This  can  be  done  by 
excluding  from  the  integrations  in  (5.1)  and  (5.o)  all  regions  of  the  (p,q) 
plane  such  that  k  »  ap  or  k  »  aq.  The  result  is  to  eliminate  spurious 
r<elaxatian  effects  on  triple  correlations  but  to  retain  in  the  energy -balance 
equation  the  transfer  effects  associated  with  strainine  of  snail  scales  by 
large  scales.  Detailed  energy  conservation  is  exactly  preserved. 
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FIGURE  CAPTION 

FIG.  1.  One-<li«ensioMi  energy  soectr\»  at  high  wavenumbers  accortling  to 
the  direct -interact ion  approxieation.  Curve  1,  unmodified  Navier-Stokes 
equation  (a  =  -).  Rj^  820.  Curve  2.  modified  Navier-Stokes  equation, 
a  =  4^,  Rj^  bOO.  Curve  3,  modifieo  equation,  o  =  2^,  R^^  ^  610.  The 
dashed  lines  show  slopes  of  -3/2  and  -6/3  for  coroarison.  The  circles 
represent  data  points  -neasurec  by  •<.  Gibson  on  the  axis  of  an  air  jet 

at  800 . 
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